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Abstract 

Driven vortex lattices have been studied in a material with strong pinning, 
such as Nb films. Samples in which natural random pinning coexists with arti- 
ficial ordered arrays of defects (submicrometric Ni dots) have been fabricated 
with different geometries (square, triangular and rectangular). Three differ- 
ent dynamic regimes are found: for low vortex velocities, there is a plastic 
flow regime in which random defects frustrate the effect of the ordered array; 
then, for vortex velocities in the range 1-100 m/s, there is a sudden increase 
in the interaction between the vortex lattice and the ordered dot array, inde- 
pendent on the geometry. This effect is associated to the onset of quasi long 
range order in the vortex lattice leading to an increase in the overlap between 
the vortex lattice and the magnetic dots array. Finally, at larger velocities 
the ordered array-vortex lattice interaction is suppressed again, in agreement 
with the behavior found in numerical simulations. 
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The influence of both apphed magnetic fleld and driving current on the type II super- 
conductors behavior is of great scientiflc and technological interest. When the vortex lattice 
is driven by an applied current in the presence of defects, it can adopt several dynamic flow 
regimes depending on the driving current and therefore on the vortex velocity. Transport 
measurements 0], neutron scattering and Bitter decoration experiments have pro- 
vided strong evidence of these different dynamics phases, including creep, plastic flow and 
ordered (elastic) flow, and even of the existence of several phase transitions as a function of 
vortex velocity . Due to this complex behavior, the vortex lattice emerges as an excellent 
model system for many interesting flelds in condensed matter physics such as self-organized 
criticality, atomic scale friction and phase transitions in the presence of disorder/order 

Theoretical works on driven vortex lattices, in the presence of either random |^ or ordered 
IO| , p!l|] pinning centers, have predicted rich phase diagrams as a function of vortex velocity 



and pinning strength. In general, most of the experimental studies have focused on dynamic 
phases in the presence of weak random pinning centers in materials such as 2if-NbSe2 single 
crystals [1^ or amorphous NbsGe fllms [jl3[ where the low values of the critical current 



allow experimental access to a wide portion of the relevant force- velocity space. However, 
driven vortex lattices materials with strong pinning, such as Nb, have received much less 
attention. In these systems, the fabrication of ordered arrays of artiflcial pinning centers by 
submicrometric lithographic techniques have proved to be a very useful tool in the study of 



vortex dynamics under tailored nanostructured pinning potentials []14-18 



In particular, depending on the size and material of the artiflcial pinning centers, the 
relative strength of the ordered and disordered pinning potentials in the sample can be tuned 
in a controlled way . 

In this work, we present results obtained from transport measurements on Nb fllms with 
ordered arrays of magnetic dots in order to study the dynamic interactions of the vortex 
lattice with random and ordered defects. Deformations in the vortex lattice frustrate the 
interactions with the artiflcial dot array for slow vortex motion, until the onset of ordering 
in the vortex lattice takes place for vortex velocities in the 1-10 m/s range. The evidence 
of this transition has been found to be independent on array geometry, indicating that the 
ordered array of dots is acting as a probe of the intrinsic behavior of the driven vortex 
lattice. Finally, for velocities above 100 m/s the effect of the ordered array on the moving 
vortex lattice disappears, in good agreement with numerical simulations. 

Magnetic (Ni) dots approximately 200 nm in diameter and 40 nm in thickness were 
prepared on Si (100) substrates using electron beam lithography combined with a lift-off 
process described elsewhere |lT5|pO|| . After this, a 100 nm thick Nb fllm was sputtered on 
top. Finally optical lithography and reactive ion etching were used to deflne a 80 fim 
wide bridge for transport measurements. Three different array geometries were fabricated: 
0.45 X 0.45yum^ square array, 0.35 x 0.5/xm^ rectangular array and a triangular array with a 
lattice constant of 0.4 /xm. 

In the case of the rectangular array, transport measurements were carried out by applying 
current parallel to the shortest lattice dimension. The dc measurements were performed in 
a hehum cryostat with a 9T superconducting magnet. The magnetic fleld is always applied 
along the fllm normal. 

As a flrst approach to the dynamics of driven vortex lattices we can consider that, in the 
presence of a driving current, the vortex lattice is set into motion with average velocity v 
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by the Lorentz force Fl = $oJ x n, where J is the transport current, $o is the quantum of 
flux, and n is the unit vector along the field direction. The vortex velocity v is determined 
by the balance between the viscous, pinning and Lorentz forces as 

r/v + av X n = Fl + (Fp) (1) 

where rj is the Bardeen-Stephen viscosity, a is the viscous Hall coefficient, and (Fp) is the 
pinning force averaged over vortex positions in time. In the dirty limit, a << t] and the 
viscous Hall force can be neglected. The resulting force-velocity curves {Fl vs v) can be 
simply derived from experimental I-V characteristics calculating the Lorentz force and using 
the Josephson relation for the electric field E = B x v, as shown in Fig. 1 for a Nb film 
with and without periodic pinning centers. For low velocities the curves approach a constant 
force value that corresponds to the static pinning force FpQ, related to the critical current 
Jc by Fpo = ^oJc- The free flux flow behavior, characterized by a linear dependence of Fl 
on V should be recovered for very high velocities, where the interaction between the pinning 
centers and the vortex lattice becomes negligible. However pure free flux flow is usually only 
accessible experimentally by fast ramp current methods, even in relatively weak pinning 
materials [Q. Often, in a wide range of vortex velocities such as plotted in Fig.|l|, the 
force velocity dependence is found to be much weaker than linear, revealing that the term 
corresponding to pinning interactions (Fp) is of the same order as the Bardeen-Stephen 
viscous force. 

The calculation of (Fp) is not straightforward, since interactions among the vortices and 
with the pinning centers induce dynamic deformations in the vortex lattice that strongly 
affect this average. The orientation of (Fp) is opposite to the vortex velocity |^T|, i.e. it 



can be written as (Fp) = — 7(f )v, with 7(f) a velocity dependent coefficient and v the unit 
vector along v. At sufiicently large driving currents, the time average of the pinning force 
is related to the spatial average by (Fp)time ~ (|Fp^|)position/|v|. The equation of motion of 
the vortex lattice is given by 

7]v + 7(t;)v = Fl (2) 

Equation (Q) implies that the effect of pinning centers can be assimilated to a kinetic viscous 
force, in a similar way to the problem of friction between two solids ^ (in this case, the 
vortex lattice and the pinning array). 

The distribution and strength of the pinning centers affect both the magnitude and 
velocity dependence of (Fp). This is clearly observed, for example, in differences between 
the two Fl — V curves in Fig.|I], that were measured in two adjacent regions of the same 
Nb film, with and without periodic pinning centers. Further insight in the behavior of the 
vortex lattice can be obtained by focusing on the differences in (Fp) depending on the 

commensurability between the vortex lattice and periodic array of dots. 

A general view of the effect of commensurability is shown in Fig.^, where we have plotted 
the voltage versus magnetic field curves obtained for different current densities in the range 
10^ - 10^ A/m^ for the Nb film with the square array of magnetic dots. Magnetoresistence 
minima are observed at equal magnetic field intervals, AHm=98 Oe. These can be explained 
by a geometric matching between the vortex lattice and the magnetic dot array. It should 
be noted that the current dependence of the magnetoresistance is not the same whether the 
measuring field corresponds to a matching condition or not. These differences show up in 
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the fact that the magnetoresistance minima are much deeper in an intermediate current 
range than at low and high currents. Actually, the critical current, shown in the inset of 
Fig.^], presents a monotonous decrease as a function of field, without any special feature 
at the matching field values. These magnetoresistance and critical current data are clearly 
indicating the existence of several different regimes: first, at f = 0, the static pinning force 
FpQ is dominated by the random pinning centers, that frustrate the effect of the ordered 
array of dots; then, in the driven vortex lattice the kinetic viscous force due to the ordered 
array of pinning centers (Fp) is enhanced in a certain velocity range, and disappears at 
large driving currents. This behavior is characteristic of small dot diameter arrays where 
periodic pinning is weaker than random pinning, and it is also observed for the rectangular 
and triangular arrays used in this study, but with Ai^m = 120 and 145 Oe, respectively, in 
good agreement with theoretical predictions of the matching conditions. 

Figures 0(a)- (c) show the force- velocity curves in the driving current range where 
commesurability minima are observed, derived from the I-V characteristics, for the Nb films 
with square, triangular and rectangular arrays of dots measured at two different fields: at 
the first matching field Hcomm for each array (filled symbols) and at a smaller field value 
Hincomm close to Hcomm but far euough to be out of matching conditions (hollow symbols). 
The obtained curves present quite similar values for low and high vortex velocities. At in- 
termediate vortex velocities (1-100 m/s) the matching curve differs from that obtained at no 
matching fields. In all cases, the overall dependence is clearly weaker than linear, indicating 
that the vortex lattice is far from the free fiux fiow regime. Two kinds of force contributions 
can be considered in this system; the sum of the Bardeen-Stephen viscous force and the 
interaction with random defects {rjv + {Fp^random)) , and the average interaction force from 
the periodic array (Fp^ordered) , which is a strongly peaked function of the magnetic field at 
the matching condition. In a first approximation, taking into account that in these sam- 
ples random pinning is much stronger than periodic pinning, {Fp^ordered) can be assumed to 
be negligible for the incommensurate vortex lattice; then the difference between the filled 
and hollow symbols curves can be used to extract the value of (Fp^ordered) at the matching 
conditions. 

This force increment, AF = F{Hcomm)-F {Hincomm) = {Fp^ordered) is shown in Figs. 3(d)- 
(f) for the three arrays considered in this work. There is a clear onset of the interactions 
with the ordered pinning array in the range 1-10 m/s, observable in the three cases. This 
onset is almost independent on the ordered array geometry, so that it must be related with 
an intrinsic change in the properties of the driven vortex lattice. The increase in (Fp^ordered) 
can be directly correlated with an enhancement in the long range order in the vortex lattice, 
since this force should be zero for a completely disordered vortex state. Therefore, the data 
in Fig. 3 show the crossover between two different regimes: first, for slow vortex motion 
{Fp^ordered) = 0, implying a plastic motion regime where the effect of the ordered array is 
frustrated by random pinning and, then, for f ~ 10 m/s the lattice flows in a more ordered 
state in which it has a better overlap with the magnetic dot array. 

In general, theoretical simulations of driven vortex lattices have predicted ||^,T^] two kinds 



of ordered states: a smectic glass phase, with quasi long range order only in the direction 
transverse to the motion, and a Bragg glass phase at higher velocities that is also ordered 
along the longitudinal direction. It is worth to note that the ordered pinning effect observed 
here occurs at Hcomm corresponding to matching between the vortex lattice with the two 
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dimensional dot array cell (this is particularly evident in the case of the rectangular array, 
where Hcomm = 120 Oe corresponds to a vortex density 5.8 x 10^ cm~^ in good agreement 
with the density of pinning centers (0.35 x 0.5/im^)~^= 5.7 xlO^ cm~^). This implies that 
both the longitudinal and transverse correlation lengths in the vortex lattice are at least as 
large as the array cell dimension when the peak in AF appears, suggesting the presence of a 
Bragg glass phase. The ordering velocity that we have found in these Nb films is two orders 
of magnitude higher than those obtained previously for other superconducting systems 



such as 2if-NbSe2 (typically of the order of 0.1 m/s). These differences can be correlated 
with the higher values of the critical current in our samples that imply stronger random 
pinning potentials [Jq ~ 10^ A/m^ in the Nb films presented in this work, while Jc ~ 10^ 
A/m^ in ^). 

Finally, it has to be noted that for vortex velocities above 100 m/s, AF decreases again, 
i.e. the interaction between the ordered array and the vortex lattice is suppressed. We 
have tried to understand this large velocity regime by modelling the moving vortex lattice 
by a ID array of point particles which move under a driving force in a combination of a 
periodic potential and random pinning centers. A similar setup, in 2D, was considered in [Q, 
although the issue of commensurability was not addressed. Our model is an approximation 



to the ID chains with short range order which are expected in a driven smectic |EBI]. The 



results are shown in Fig.^. Commensurability effects appear only at one matching field, due 
probably to the length of the chain used (20-30 vortices). It is interesting to note that, as 
in the experimental data shown in Fig. 2, these effects dissappear at large driving currents. 
The simulations show that there are, at least, three different regimes: i) At low driving 
currents, the lattice remains disordered at all fields, leading to structureless V-B curves, ii) 
At intermediate currents, the vortex lattice shows some degree of order, which increases as 
a function of the field. In this regime the minimum in the V-B curves at the matching field 
is most pronounced, iii) Finally, at high driving currents, the lattice is very well ordered, 
with no substantial dependence on the applied field. In this regime, the effects at Hcomm are 
very weak, if any. It is worth to note that the same effect of reduction in the interaction 
force at high velocities has been numerically predicted in simulations of the similar problem 
of friction between two lattices at atomic scale [Q. 

In summary, we have analyzed the force-velocity characteristics for several Nb films with 
square, triangular and rectangular arrays of Ni dots, in which ordered and disordered defects 
coexist. By comparison of the curves at matching magnetic field and away from matching 
conditions we have found that there is a clear increase in the interactions between the vortex 
lattice and the ordered array of dots in the vortex velocity range 1-100 m/s. The sudden 
enhancement in (Fp^ordered) has been attributed to the onset of quasi long range order in 
the vortex lattice at this driving velocity. This result is independent on the array geometry, 
indicating that a change in the intrinsic properties of the vortex lattice is being observed. 
For very high vortex velocities, commensurability effects are suppresed again, in agreement 
with numerical simulations. 
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FIG. 1. Force vs velocity curves derived from the I-V characteristics for (a) a simple 
Nb film and (b) a Nb film on a rectangular array of Ni dots, grown on the same substrate, 
at H= 81 Oe and T = 0.98T(7. This field corresponds to the first matching field of the 
rectangular array. 

FIG. 2. (a) V — B curves obtained on the square array at T = 0.99Tc for different current 
densities (J = 6.25 x 10^1.25 x 10^3.1 x 10^ 6.25 x 10^9.4 x 10^ 1.25 x 10^ 1.56 x 10^ 
and 1.87 x \Q^A/w? for the (a)-(h) curves respectively). Inset shows the field dependence 
of the critical current at the same temperature. 

FIG. 3. (a) Force-velocity characteristics for (a) a Nb film on a 0.45 x 0.45/im^ square 
array of dots (filled symbols, Hcomm = 98 Oe; hollow symbols, Hincomm = 60 Oe) (b) a Nb 
film on a triangular array of dots of lattice constant 0.4/xm (filled symbols, Hcomm — 145 Oe; 
hollow symbols, Hincomm = 114 Oe) (c) a Nb film on a 0.35 x 0.5/im^ rectangular array of 
dots (filled symbols, Hcomm — 120 Oe; hollow symbols, Hincomm — 90 Oe); (d), (e), (f) Force 
increment AF = F{Hcomm) — F{Hincomm) (i-e. from no matching conditions to matching 
conditions) as a function of vortex velocity for the Nb films in (a), (b) and (c) respectively. 
Lines are guides to the eye. 

FIG. 4. V — B characteristics of a model of moving vortices under the infiuence of a 
periodic potential and random pinning (see text), corresponding to: filled symbols, J = 10 Jo 
and hollow symbols, J = 2.5 Jo, where Jo is the critical current at 0.6 Hcomm- 
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FIG. 1. M. Velez et al. 
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FIG. 2. M. Velez et al. 
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FIG. 3. M. Velez et al. 
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FIG. 4. M. Velez et al. 
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